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1. INTRODUCTION 

Many advances of chaos theory and chaotic systems have been actively carried out in the last few 
decades [1]-[3]. Classical examples of 3-D chaotic systems include the Rikitake dynamo system [4], and Liu 
system [5]. New chaotic systems have been also reported in the literature such as Sambas systems [6]-[9], 
Chen system [10], and Sprott system [11]. Chaos theory has many applications of nonlinear chaotic 
oscillators [12]-[15]. Idowu et al. [13] discussed the adaptive control and circuit implementation of a new 3- 
dimensional chaotic system with quadratic, cubic and quartic nonlinearities. Hu et al. (2 discussed adaptive 
control methods for Chua’s chaotic system [15]. 

Chaotic systems have applications in neuron models [16]-[18]. Akaishi et al. [16] presented a new 
theoretical model for the nonlinear clinical and pathological manifestations in multiple sclerosis. 
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Luo et al. [17] presented new results for the adaptive passive control of the FitzHugh-Nagumo chaotic 
neuron model. Hong [18] used adaptive control to discuss the chaos in neurons and devise synchronization 
schems for Hindmarsh—Rose neuron model. 

Chaotic systems have applications in neural networks [19]-[21]. Akhmet and Fen [19] discussed the 
generation of cyclic and toroidal chaos by Hopfield neural networks. Wang et al. [20] derived new results for 
the synchronization of fuzzy cellular neural network attractors via adaptive control method. Bao et al. [21] 
derived new results for memristive neural networks with threshold electromagnetic induction. 

In order to enable synchronisation of chaotic systems for encryption schemes and secure 
communication devices, many useful control techniques have been devised in the literature such as adaptive 
control [22], [23], active control [24], [25], fuzzy control [26], [27], backstepping control [28], [29], and 
sliding mode control [30], [31]. Sliding mode control (SMC) is a special robust approach for variable 
structure control (VSC) knowledge. Since SMC technique-based control law is effective and guarantees both 
the occurrence of sliding motion and controlling of the nonlinear systems [32], [33]. So, in this paper, an 
improved SMC approach is proposed for the synchronous control of chaotic system with Lyapunov-based 
SMC strategy. The chattering problem, which appeared in these designed sliding mode control laws, is 
handled thanks to a proposed quasi sliding mode control (QSMC) technique [34], [35]. 

In this research work, quasi-sliding mode control (QSMC) approach to design the synchronous 
control of chaotic system with Lyapunov-based SMC strategy is proposed. In this work, a QSMC based 
master and slave Rikitake chaotic circuits are proposed to guarantee asymptotic synchronization. QSMC 
design has advantages of fast convergence, robustness and unperturbed by system uncertainties. MATLAB 
plots are exhibited to illustrate the efficiency of the proposed QSMC control scheme. 


2. RIKITAKE CHAOTIC DYNAMO SYSTEM 

Explaining research chronological, including research design, research procedure (in the form of 
algorithms, Pseudocode or other), how to test and data acquisition [1]-[3]. The description of the course of 
research should be supported references, so the explanation can be accepted scientifically [2], [4]. Rikitake 
chaotic dynamo system (1958) is described as follows [4]: 


x= —bx + zy 
ý = —by + (z —- a)x (1) 
ż = 1- xy 


We use the notation X = (x, y,Z) to represent the state of the Rikitake system (1). Here, a,b are positive 


parameters. 

The Rikitake dynamo system is a famous model that is capable of executing nonlinear chaotic 
oscillations similar to the chaotic oscillations as revealed by palaeomagnetic study. The behaviour of two 
disk dynamos coupled to one another was examined by Rikitake [4] in relation to the earth's magnetic field. It 
was found by Rikitake [4] that reversals of electric current and magnetic field occur in the Rikitake two-disk 
dynamo system (1) unlike the case of a single disk dynamo. 

In this work, we assume that (a,b) = (3,1) and X(0) = (1,1,1). The Lyapunov exponents of the 
Rikitake system (1) are obtained using MATLAB as: 


L, = 0.1661, L, = 0, L = —2.1661 (2) 


Since L, > Oand L, +L, + L} < 0, we conclude that the Rkikitake system (1) is chaotic and dissipative. 
Thus, the system orbits of the Rikitake dynamo system (1) are ultimately confined into a specific limit set of 
zero volume and the asymptotic motion settles onto a chaotic attractor. The Kaplan-Yorke dimension of the 
Rikitake dynamo system (1) is found as; 


Ly4+L2 
|L3| 





This shows the high complexity of the Rikitake dynamo system (1). We find that the Rikitake system (1) is 
invariant under the coordinates transformation, 


(x,y,z) > (=x, =y, 2) (4) 


The invariance holds for all values of the parameters aand b. This shows that the Rikitake dynamo 
system (1) has rotation symmetry about z-axis. The equilibrium points of the Rikitake dynamo system (1) are 
easily obtained as: 
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1.8174 —1.8174 
E =|0.5503|, Ep = -0.503| (5) 
3.3028 3.3028 


A simple calculation shows that F, and F,are saddle-foci equilibrium points. Hence both are 
unstable. Thus, the Rikitake dynamo system exhibits self-excited chaotic oscillations. The phase portraits of 
the Rikitake dynamo system (1) are showed in Figure 1. From this figure, we see that the Rikitake dynamo 
system (1) exhibits a double-scroll chaotic attractor. 














(d) 


Figure 1. Numerical simulations of the Rikitake dynamo system (1) for (x(0), y(0), z(0)) = (1,1,1) and 
(a,b) = (3,1): (a) x - y plane, (b) y - z plane, (c) x - z plane, and (d) R? 


3. CIRCUIT IMPLEMENTATION OF THE RIKITAKE DYNAMO SYSTEM 

Circuit design of chaotic systems has important applications in engineering [36], [37]. Zhang and 
Liao [38] discussed chaos synchronization and chaos in coupled memristor-based FitzHugh-Nagumo circuits 
with memristor synapse and they gave a design of the memristor-based circuit of FitzHugh-Nagumo model. 
The Rikitake dynamo system (1) is designed in Multisim by making use of the electronic components such as 
capacitors, resistors, operational amplifiers TLO82CD and analog multipliers AD633JN. 

For the Multisim circuit design, operating voltage of the operational amplifier is +15V. The x,y,z 
signals of the Rikitake dynamo system (1) are emulated by the output voltages of the integrators. For the 
circuit implementation of the Rikitake system (1), we adopt the following scaling of the signals x, y, z of the 
Rikitake system (1): 


X= Z2 
Y =2y (6) 
Z=22 


TELKOMNIKA Telecommun Comput El Control, Vol. 19, No. 4, August 2021: 1428 - 1438 


TELKOMNIKA Telecommun Comput El Control 


The rescaled Rikitake dynamo system is given by; 
X =-bX+5YZ 


Ý = -bY + (+Z - a) x 
Z=2-=XY 


We apply the Kirchchoff’s laws to derive the circuit presented in Figure 2 as follows; 


1 1 

















X =-—X+ YZ 
Ci C1R> 
Y = -—_Y + xz -—_ x 
C2R3 C2R4 C2R5 
Fe aay 
C3R7 C3R6 


We take values of the system components as follows: 
Let C, = Cz = Cg = 1nF. 

Let R, = R3 = 400k andR , = R; = Rs = 800K. 
Let Rs = 133.33 kQand R, = 200 kf. 

Let Rg = Ro = Rip = R41 = 100 kR. 


o 


(7) 


(8) 
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The values of various capacitors and resistances were assumed using parameter values of the 
Rikitake dynamo system and practical considerations. We obtain Multisim outputs of the Rikitake circuit 
which are shown in Figure 2. It is clear that the Multisim outputs in Figure 3 exhibit the double scroll 
attractor of the system (7). Thus, the system (7) undergoes chaotic behavior. The Multisim outputs of the 
Rikitake circuit (7) in Figure 3 are consistent with the MATLAB plots given in Figure 1 of the Rikitake 


dynamo system (1). 


Figure 4 illustrates the Fourier spectral analysis plot in Multisim for the chaotic signal X of the 
Rikitake circuit (7). The frequency range is 5 kHz, maximum peak 750 Hz. It corresponds to a prevailing 
frequency of the implementing oscillating loop. The power spectra of the produced signals are broadband, 
typical of chaotic signals. The frequency is quite low. Thus, the new chaotic system can only be used for low 
frequency applications. Fourier spectral analysis plot of Figure 4 confirms chaotic oscillations in the Rikitake 


dynamo system (1). 
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Figure 2. The electronic circuit schematic of the Rikitake dynamo circuit (7) 
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Figure 3. Multisim outputs of the scaled Rikitake dynamo circuit (7) in 
(a) x - y plane, (b) y — z plane, (c) x — z plane 
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Figure 4. Fourier spectral analysis plot in Multisim for the signal X of the Rikitake circuit (7) 


4. QUASI-SLIDING MODE CONTROL BASED SYNCHRONISATION OF RIKITAKE 
CIRCUITS 

Sliding mode control (SMC) is proposed to tackle many control problems in engineering 
applications [39]. There are many variants of SMC such as higher order SMC [39], super-twisting SMC [40], 
integral SMC [41], robust integral SMC [42], and terminal SMC [43]. In this section, quasi-sliding mode 
control (QSMC)-based master and slave Rikitake chaotic dynamo systems are proposed to guarantee 
asymptotic convergence of the synchronization errors of the master-slave Rikitake circuits. The parameters of 
the Rikitake dynamo systems are taken as in the chaos case, i.e. (a, b) = (3,1). 
Master system: 


Vin = —Ym + (Zm — 3)Xm (9) 


fiz = —Xm + ZmYm 
Zm = 1 — XmYm 


Slave system: 


Vs = =y F (25 — 3)X5 (10) 


i = —X; + Z;Y⁄s tu 
Ze Imay 


In the slave Rikitake system (10), U is the quasi-integral sliding mode control (QSMC). We next describe the 
synchronising errors between the Rikitake systems as follows; 


e1 (t) = xs (t) — Xm (t) 
e2(t) = Yst) — Ym (t) (11) 
e3 (t) = Zs (t) — Zm (t) 
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The error dynamic system is calculated as follows; 


é (t) = —e1 + ZsYs — ZmVm + U 
€2(t) = —e2 + (Zs — 3)Xs — (Zm — 3) Xm (12) 
€3(t) = —XsVs + XmYm 

The proposed quasi-sliding mode control (QSMC) U can be defined by; 


S 





E (13) 
The switching surface can be chosen as follows: 

S = e, + À& (14) 
Here, A, ô, ware constants and 

n = |[-e2 + (Zs — 3) Xs — (Zm — 3)Xm] + Al—e1 + ZsYs — ZmYmll| (15) 


The quasi-sliding mode control (QSMC) U is designed to derive the error dynamics (11) satisfying 
the reaching condition s(t)s(t) <0 for the sliding surface s(t) =O. The invariance conditions of the 


sliding manifold are given as follows: 

s(t) = 0 and s(t) = 0 (16) 
To simplify the notation, we define 

A = —é€, + Z;Ys — ZmYmand B = —e, + (Z; — 3)x; — (Zm — 3) Xm (17) 
In (15) can be then simplified as follows: 

n=|B+AA| (18) 
We also note that; 

S=é@,+Ae,; = B+A(A +u) (19) 
That is, 

$=B+AA+Au (20) 
The proposed Lyapunov function can be provided as; 


V ==s? (21) 


We calculate the time-derivative of V as follows: 


V=ss=s(B+AA+Au) = s(B + AA) + sdu (22) 


Thus, it follows that 











. wns? _ _ = |s|6 
V <n|s|— =" = nls| — wn (Isl - 35) (23) 
We note that 
|s|o 
|s|+5 — (24) 


Thus, we can simplify (23) as follows: 


V < nls] — wn(|s| - 6 = (1 -= wn (Isi - #5) (25) 


We define 
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dg = wet (26) 


With w > land | s(¢) |> p, we can conclude that 


V =(1-—w)(|s| — 5g) < 0 (27) 


By Lyapunov stability theory [44], we conclude that s(t) — Oasymptotically as to. By 
quasi-sliding mode control theory, we have established the chaos synchronization between the master 
Rikitake system (8) and slave Rikitake system (9). 


4.1. Numerical simulations 

In this section, we give MATLAB simulation results for the asymptotic synchronization between the 
master Rikitake system (8) and slave Rikitake system (9). The QSMC control parameters are set as A = 30, 
ô = 0.003and w = 1.5. The initial data for the master system (8) is taken as x,,(0) = 0.3, y,,(0) = 0.2 and 
Zm(0) = 0.1. The initial data for the slave system (9) is taken as x,(0) = 0.6,y,(0) = —0.1 and 
Z,(0) = —0.4. The synchronization of the respective state trajectories between the master and slave Rikitake 
dynamo systems (8) and (9) is shown in Figure 5. The synchronization error vector 
e(t) = (e (t), e2(t), e3(t)) between the master and slave Rikitake circuits is depicted in Figure 6. 
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15 
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Figure 5. Complete synchronization of the master and slave Rikitake systems (8) and (9) 
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Figure 6. Time-plot of the synchronization error e(t) = (e1 (t), e2(t), e3(t)) between the master and slave 
Rikitake dynamo systems (8) and (9) 
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5. CONCLUSION 

In this paper, we gave a brief study of the dynamic properties of the Rikitake circuit (1958) and 
presented a circuit design for the Rikitake chaotic attractor using Multisim. The main result of this paper is 
the chaotic synchronisation of Rikitake circuits using a quasi-sliding mode control (QSMC) design. Suitable 
quasi-sliding mode control controllers were designed by using Lyapunov stability theory. We designed the 
switching surface for the synchronisation error dynamics. Then a quasi-sliding mode control was considered 
to guarantee the synchronisation effect on the master and slave Rikitake circuits. MATLAB simulations were 
shown to illustrate the signal plots of Rikitake circuits and synchronisation results via QSMC. 
Synchronisation of chaotic systems and circuit design have several applications in engineering areas such as 
secure communication systems and encryption devices. 
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